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ABSTRACT 



To tackle the sign problem in the simulations of systems having indefinite or 
complex- valued measures, we propose a new approach which yields statistical errors 
smaller than the crude Monte Carlo using absolute values of the original measures. 
The ID complex-coupling Ising model is employed as an illustration. 

1. The sign problem 

If the measure p(x) of a generating function suffers from the sign fluctuation 
then another positive definite function p(x) must be chosen for the Monte Carlo 
(MC) evaluation of the expectation value of an obsevable. With the choice 



the sign of p(x) is treated as part of the quantity whose expectation is to be mea- 
sured; hence the name average sign approach 1 . However, the fluctuation of sign 
of the measure over configuration space renders ineffective the sampling guided by 
this crude MC method; this is the content of the sign problem. Many approaches 
have been proposed to tackle the sign problem but none is satisfactory 2 . 

2. Our new approach 

We can write the generating function as integrals over two configurational sub- 
spaces, J x p(x) = J x J Y p(X } Y) } in such a way that the multi-dimensional integral 
over Y can be evaluated analytically or well approximated: 







One can easily prove by variational techniques that the weight 
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is the MC weight that minimises the variance 
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Z being the complex partition function. It then follows that the variance for this 
new weight is not bigger than that for p(X,Y). 
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the second line follows from (1) and (3); the last line is from (2) and always less 
than or equal to zero because of the triangle inequality. The equality occurs if and 
only if p(X } Y) is semi-definite (either positive or negative) for all X-configuration. 
In particular, when there is no sign problem in the first place, expression (4) yields 
the same statistical deviation as the crude one. 

To deal with complex integrands, of which indefinite measures are special cases, 
we adopt the definition above of the variance extended to the absolute values of 
complex numbers. Statistical analysis from this definition is the same as the stan- 
dard analysis; except that the range of uncertainty should now be depicted as the 
radius of an 'uncertainty circle' centred on some central value in the complex plane. 
All the manipulations above also remain valid. 

Our approach is now clear. The measures are first summed over a certain sub- 
space, the integration (2) above, to facilitate some partial phase cancellation. Ab- 
solute values of these sums are then employed as the MC sampling weights (3). 



3. Complex coupling Ising spins 



Owing to the nearest-neighbour interactions, the lattice can be partitioned into 
odd and even sublattices, of which the Ising spins s 8 - (= ±1) on site z(G the sublat- 
tice) do not interact with each other. Absolute values of sums of the complex- valued 
weights over the even sublattice, say, are the new Monte Carlo weights, with complex 
coupling J and no external field 

Q(M) ~ n |cosh(J(s,- + s,- +1 ))|. (6) 

odd sites 

In our simulations, periodic boundary condition is imposed on the ID chains 
of 128 Ising spins which become 64 spins after the partial summation. Ensem- 
ble averages are taken over 1000 configurations, out of 2 128 possible configurations. 
They are separated by 30 heat bath sweeps which is sufficient for thermalisation 
and decorrelation in all cases except perhaps one, as will be demonstrated shortly. 
A heat bath sweep is defined to be one run over the chain, covering each spin in 
turn. The numbers of trials per sweep are different before and after the partial 
summation because the numbers of spins to be updated are not the same. Both 
hot and cold initial configurations are used. We present in the table some measure- 
ments by exact, improved MC and crude MC methods respectively. Expressions 



for the magnetisation and susceptibility are obtained from appropriate derivatives 
of corresponding partition functions. The autocorrelation of the unit operator at 



Coupling 


Magnetisation per spin 


Susceptibility per spin 


(0.1,0.1) exact 
improved 
crude 


(0,0) 

(0.0007,0.0000) [0.0024] 
(-0.0032,0.0016) [0.0062] 


(1.1971,0.2427) 
(1.1921,0.2459) [0.0356] 
(1.1651,0.2883) [0.1460] 


(0.01,0.1) exact 
improved 
crude 


(0,0) 

(-0.0006,-0.0001) [0.0021] 
(0.0023,-0.0080) [0.0055] 


(0.9999,0.2027) 
(1.0481,0.2231) [0.0268] 
(1.0868,0.1374) [0.1198] 


(0.01,0.5) exact 
improved 
crude 


(0,0) 

(-0.0031,-0.0047) [0.0027] 
(-0.0522,0.0820) [0.2074] 
(0.0123,-0.0723) [0.0830] 


(0.5512,0.8585) 
(0.5604,0.8163) [0.0428] 
(0.6808,5.4507) [9.2788] 
(1.5323,1.3885) [2.3122] 



J = (0.01,0.5) show that the noise is too overwhelming in the crude simulation to 
tell whether 30 sweeps are sufficient for thermalisation or not, hence both the hot 
and cold starts are shown in the table. In contrast, the improved simulation is very 
well-behaved. 



4. Conclusion 

Our approach can offer substantial improvements over the crude average sign 
method and may work even when the later fails, in the region of long correlation 
length and vanishing partition function. The choice for splitting of the integration 
domain is arbitrary and its effectiveness depends on the physics of the problem. If 
the interactions are short-range (not necessary nearest-neighbour), maximal, non- 
interacting sublattices can always be chosen to provide a natural splitting, which is 
the particular splitting for our illustrative example. 
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